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Rheological properties, especially 'shear-thinning', of dense colloidal dispersions are discussed on 
three different levels. A generalized phenomonological Maxwell model gives a broad framework con- 
necting glassy dynamics to the linear and non-linear rheology of dense amorphous particle solutions. 
First principles mode coupling theory calculations for the time or frequency dependent shear modu- 
lus give quantitative results for dispersions of hard colloidal spheres in the linear regime. Schematic 
models extending mode coupling theory to the non-linear regime recover the phenomenology of the 
generalized Maxwell model, and predict universal features of flow curves, stress versus shear-rate. 
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I. INTRODUCTION 



The rheological behavior of dispersions of rigid col- 
loidal particles is immensely rich 0, Universally 
they exhibit shear-thinning, that the viscosity r\ de- 
creases from the Newtonian-plateau upon increasing the 
shear rate 7. Another phenomenon in dense colloidal 
dispersions, or dispersions of particles with strong at- 
tractions, is that they undergo arrest-transitions and 
solidify to (mostly weak) amorphous solids (glasses) 
HH E HQ H In this contribution, we want to 
establish under which circumstances and simplifications 
one may expect that the physics of the glass transition 
determines or even dominates the rheological behavior 
of colloidal dispersions. Our discussion supplements a 
recent comparison of theoretical computations and mea- 
surements of flow curves in colloidal model systems where 



this connection is explored quantitatively 



syste 
10]. 



While recently microscopic approaches to shear- 
thinning in colloidal dispersions have been presented 
[Til ll2L Il3l Pill ], in this contribution we first argue on 
a more qualitative level and sketch a phenomenological 
picture, that goes back to Maxwell. Then, we discuss in 
how far simple 'schematic models' derived from the first 
principles approaches in Refs. 0,0,0,0, support 
the generalized phenomenological Maxwell model. The 
concepts we discuss build upon insights into the struc- 
tural dynamics obtained with the mode coupling theory 
(MCT) of Gotze and coworkers [l5|, Hg, an d thus ex- 
tend a microscopic approach which has amply been tested 
and compared to computer simulation studies and exper- 
imental data from quiescent model systems . We also 
present MCT predictions for the linear viscoelasticity of 
colloidal hard spheres which may be tested in future sim- 
ulational or rheological studies. 



II. A PHENOMENOLOGICAL GENERALIZED 
MAXWELL MODEL 

A. Quiescent system 

Before discussing the connection of glassy dynamics 
with the rheology of dense colloidal dispersions, it ap- 
pears helpful to reconsider the oldest phenomenological 
description of the glass transition which goes back to 
Maxwell. Maxwell considered the most general consti- 
tutive equation relating (transversal) stress a and shear 
rate 7 within the framework of linear response theory: 



a(t) 



dt' G to \t-t') j(t') . 



(1) 
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Here, G tot (t) is the shear modulus and is, within linear 
response, independent of the applied shear rate 7, which 
in the most general situation is a time-dependent func- 
tion. This relation is exact for the linear response of 
systems with free energy corresponding to a fluid phase 
provided that the system was in equilibrium before the 
(smoothly switched-on) application of shear 0. For 
colloidal solutions, Batchelor was the first to discuss the 
three contributions to the stress which arise: in the sol- 
vent, from the potential interactions among the particles, 
and from the solvent particle interactions |l9j. For the 
dispersions of interest here, the solvent can be consid- 
ered a continuum fluid characterized by a know solvent 
viscosity rj s , so that G solv (i) = rj g S(t). 

Solvent-particle interactions (viz. the HI) act instan- 
taneously if the particle microstructure differs from the 
equilibrium one, but do not themselves determine the 
equilibrium structure. If one assumes that glassy arrest 
is connected with the ability of the system to explore 
its configuration space and to approach its equilibrium 
structure, then it appears natural to assume that the 
solvent particle interactions are characterized by a finite 
time scale thi- In a most naive way, the contribution 
from HI, G m (t), to the shear modulus G tot (t) would thus 
be written as: G m (t) = G™ exp {— t/mi}, where the 
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time thi stays finite and varies little with the distance to 
glassy arrest. Considering time-independent shear rates, 
7(t) = 7, so that 

/>oo 

,7 = 777 = 7 / dtG tot (t), (2) 
Jo 

HI would thus lead to an increase of the high frequency 
viscosity above the solvent value; this value shall be de- 
noted as = rj s + G™ r H i. 

Maxwell in his studies of viscous liquids considered 
the potential contribution, to be abbreviated as G(t), 
to the shear modulus, which arises from the forces be- 
tween the particles. Its microscopic definition is, G(t) = 
(o'xyi^^xyiO)) /(ksTV), where the relevant element of 
the fluctuating stress tensor is given with the poten- 
tial forces Fj and particle coordinates by cr xy (t) = 
— J2 i Ff (t) yi (t) . Kinetic contributions that would dif- 
fer between Newtonian (which Maxwell considered) and 
Brownian (the colloidal dispersions considered here) sys- 
tems do not contribute for symmetry reasons. He pos- 
tulated that close to the glass transition the modulus 
develops a slow relaxational process that is characterized 
by a time scale r, which grows beyond bounds: 

G(t) =Goo exp{--} . (3) 

T 

The amplitude of the slow relaxational process, Goo, can 
be considered an almost constant that changes much 
more slowly with density or other parameters than the 
time scale r. The divergence of r upon solidification 
often is modeled by a Krieger-Dougherty, t/tq cx 
(1 — <fi/<f) m )~ 2 , or a Doolitlc form, lnr/ro oc (1 — (/>/0 m ) -1 . 
Here the packing fraction at the maximal fluid density 
is denoted as </> m , and the time-unit tq should be con- 
nected to the viscosity 7700 that a particle feels because of 
the surrounding solvent. The slow process described by 
Maxwell, which in the glass literature is called a-process, 
causes a marked and strongly density dependent increase 
of the total viscosity over 7700 : performing the integration 
gives 770 = t]oo + Goo t. In this picture, the total shear 
modulus as a function of time consists of the three men- 
tioned terms: G tot (t) = j] s 8{t) + G HI (i) + G(t). 

The naive picture sketched here, is not correct for a 
number of reasons. It is well known that for hard spheres 
without HI the shear modulus diverges for short times, 
G HSnoIiI (t ->• 0) ~ t- 1 / 2 . Lubrication forces, which keep 
the particles apart, eliminate this divergence and ren- 
der G tot (t -> 0) finite 20]. Thus, the simple separa- 
tion of G m (t) and the potential G(t) is not possible for 
short times, at least for particles with a hard core. More- 
over, comparison of simulations without and with HI has 
shown that the increase of (770 — 7700 ) /Voo depends some- 
what on HI, and thus not just on the potential interac- 
tions as implied. 

Nevertheless the sketched picture provides the most 
basic view of a glass transition in colloidal suspensions, 
connecting it with the increase of the structural relax- 
ation time r. Increased density or interactions cause 



a slowing down of particle rearrangements which leave 
the HI relatively unaffected, as these solvent mediated 
forces act on all time scales. Potential forces dominate 
the slowest particle rearrangements because vitrification 
corresponds to the limit where they actually prevent the 
final relaxation of the microstructure. The structural re- 
laxation time diverges at the glass transition, while r HI 
stays finite. Thus close to arrest a time scale separation 
is possible, r 3> r HI . 

While the parameter Goo in Maxwell's model might 
appear to play a rather boring role as amplitude of the 
final relaxation process, its actual significance is quite 
high. It gives the elastic constant of the arrested struc- 
ture in the limit of 1/t = (viz. in the glass), or of the 
glassy structure present on intermediate times, t <C T. 
Thus in the glass or for shear rates j(t) that vary more 
rapidly with time than the modulus, a Hookian elastic 
response follows from Eq. QJ: 

a(t) = Goou(t) + dissipative terms , where (4) 
u(t) = I dt'j(t') forl/r = 0, (5) 

where u(t) is the strain imposed on the amorphous solid. 
Clearly, Goo plays the role of a spring constant for the 
elastic restoring forces in the solid which grow with the 
imposed distortion measured by the strain, viz. u(t) = 

fdt' dv x (t')/dy. 

An interesting difference from this hall mark behavior 
of glasses to colloidal coagulation or gelation at low pack- 
ing fractions should be noted. Experiments on a wide 
range of gels at high dilutions have shown that the elas- 
tic constant Goo vanishes at the colloidal gelation line; it 
then grows with a strong power-law upon entering deeper 
into the gel phase |21|. It has been pointed out that 
this behavior can be explained by an underlying percola- 
tion like transition. Thus the observed 'Trappe-scaling' 
points to a deep difference between colloidal coagulation 
at very low packing fractions and colloidal vitrification 
at higher ones. At the latter transition, Goo is finite, as 
first postulated by Maxwell. Interestingly, this difference 
does not show up in universal jamming phase diagrams. 
Intriguing data have been obtained for intermediate con- 
centration range where both pictures transform, but fur- 
ther studies appear very fruitful |2^. 

The elastic constant Goo serves to indicate a difference 
of the present Maxwell's view on the glass transition to 
another colloid theory. It is the approach by Brady [23| . 
who identifies the solidification of amorphous dispersions 
with random close packing, <f> m = (f> IC p. In this approach, 
the modulus G tot relaxes in a single relaxational process 
characterized by the time scale r s = a / a is the par- 
ticle radius, and Dq the density dependent short-time 
self diffusion constant. While Brady's final result for the 
viscosity, 770 oc (1 — <f)/(f) m )~ 2 , is similar to glassy behav- 
ior (MCT predicts a divergence 770 oc r cx (1 — <j)j '4> c )~ 2 - 6 
for hard spheres), part of the divergence of r/o comes 
from a divergence of the elastic constant upon approach- 
ing random close packing, Goo oc (1 — <p/(p m ) . The 
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time scale r s also diverges at solidification, but con- 
tributes only comparable to the amplitude of the re- 
laxational process, r s oc (1 — (fr/ifim)" 1 . Thus, a quite 
different scenario of arrest at random close packing to 
arrest at the glass transition should become apparent if 
G tot (t) or its Fourier-transforms, elastic G'(ui) and loss 
G"(u>) moduli are measured. (We use the convention 
G"\J) + iG'{uj) = uj ffite wt G{t)) While at the glass 
transition the dynamics should only become slower, at 
random close packing the amplitude of the stress at in- 
termediate times should also increase. Besides this qual- 
itative difference, a quantitative difference also arises in 
the values of arrest, namely the glass transition density, 
called <j) c in MCT and measured in, lies below random 
close packing: <j) c w 0.58 < <j) rcp w 0.64 0, S H 0, 0, @ . 
Particles in a glass thus still have local volume to re- 
arrange; they can explore their respective 'cages', but 
cannot relax the overall microstructure. 



B. Under shear 

The utility of the phenomenological Maxwell model of 
the glass transition is further stressed by the ease by 
which it can be extended to the shear-thinning of col- 
loidal fluids and yielding of colloidal glasses under shear. 
The central concept required is that shearing speeds up 
the relaxation of fluctuations and kills off the long-time 
memory present in very viscous fluids. This requires that 
shear causes the slow structural or a-relaxational pro- 
cess in G(t) of the Maxwell model to decay faster. It is 
natural to connect the relaxation time of the structural 
relaxation to the inverse of the shear rate, and thus to 
formulate on the phenomenological level: 

g^(t,j) = G tot (t) exp{-*|7|}. (6) 

Here, the generalized non-linear shear modulus g(t, 7) 
is a function describing the transient fluctuations in a 
dispersion sheared from time t = on. Equation © 
expresses that under shear transient stress fluctuations 
either decay via the relaxation channels present without 
shear, or via a shear induced relaxation mechanism whose 
intrinsic decay rate is set by the shear rate itself. Pos- 
tulating that the shear rate enters in a non-analytic way 
via the absolute value into the phenomenological descrip- 
tion is of course, a priorL questionable, and requires the 
derivation given in Refs. [Ulllj]; see also section IV. The 
result Eq. © is quite satisfying, because, intuitively, one 
expects the direction of the solvent flow profile to be ir- 
relevant for the magnitude of the stress or viscosity even 
in the non-linear regime. 

The generalization of Eq. QJ, that connects 5(^7) to 
the steady state stress in the case of constant shearing is: 

/>oo 

a(j) = 77(7)7 = 7 / dtg* M ™(t,j) 
Jo 

= 7 <Vs + T—|T| + ... > ■ 7 

I 1 + 7 Thi 1 + 7 t J 



This is the final expression of the extended Maxwell 
model for the stationary viscosity or stress as a function 
of the shear rate, 77(7). Variations of the packing frac- 
tion or of other control parameters of the system without 
shear enter via the Maxwell relaxation time r, and con- 
sequently 770, which may for example be modeled by a 
Krieger Dougherty expression below the glass transition 
density, and by r = 00 = 7/0 above it. All other constants 
appearing explicitly or implicitly in Eq. J7J may, in the 
most simple approach, be taken as control-parameter in- 
dependent constants. 

The discussed phenomenology is quite rich with flow 
curves that exhibit (i) a Newtonian viscosity plateau 
given by 770 at low shear rates and below the glass tran- 
sition (for I7I <c 1/r <C 1/thi), (m) a shear thinning 
region for 1/r <C I7I <C 1/thi below the glass transition, 
(Hi) a solid-like yield flow for = 1/r -C |7| <C 1/thi 
above the glass transition, and (iv) a second Newtonian 
plateau characterized by 77^ following the shear thinning 
region while still I7I <C 1/thi holds. In this region of 
high shear, where the HI are dominant, the suggested 
modeling of G m (t) is not based on theoretical insights 
from our microscopic approach, and thus is not reliable. 
Therefore, the most sensible approach, that we also will 
follow in the remainder of this presentation, is to restrict 
the shear rates to the region |7|rni <C 1. Then Eq. Q 

simplifies to 77(7) = t;^ + + j^r^ , which captures 
the small shear rate limit close to vitrification, where 
r 3> tbi, and 770 S> 77,30. The notation a* has been in- 
troduced, to stress the dimensionalities; from Eq. Q 
c* = Goo follows, while in more general settings a* is 
more cl osely related to the critical yield stress at solidi- 
fication [14] . 

Aim of section IV is to present and discuss a schematic 
model, which can be derived within our microscopic ap- 
proach [Til IT^ , and which recovers the phenomenology 
just described. 



III. VISCOELASTICITY OF QUIESCENT 
COLLOIDAL FLUIDS 

Mode coupling theory (MCT) appears to capture the 
cage effect and predicts that it dominates the slow relax- 
ation of structural correlations [16j . Density fluctuations 
play an important role because they are well suited to 
describe the structure of the particle system and its re- 
laxation. Moreover, stresses that decay slowly because 
of the slow particle rearrangements, MCT argues, also 
can be approximated by density fluctuations using effec- 
tive potentials. Density fluctuations not at large wave- 
lengths, but for wavelengths corresponding to the average 
particle distance turn out to be the dominant ones. In 
agreement with the picture of the caging of particles by 
structural correlations, the MCT glass transition is inde- 
pendent on whether the particles move ballistically in be- 
tween interactions with their neighbors (say collisions for 
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hard spheres) or by diffusion. Structural arrest happens 
whenever the static density correlations for wavelengths 
around the average particle distance are strong enough. 
Thus, it had been realized early on that this basic version 
of MCT, the 'idealized MCT', that was developed for a 
description of the structural arrest in atomic systems also 
needs to apply to colloidal dispersions 0] . 



T — T 
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FIG. 1: Equilibrium potential shear modulus G(t) (in units of 
fcajT/o 3 ) for colloidal hard spheres with radius a for packing 
fractions close to solidification at <f> c . Densities are measured 
by the separation parameter s = ((/> — <f> c )/<f) c = ±10~' n ' / ' 3 , 
and labels denote the value n. Positive values belong to glass 
(e > 0), negative to fluid states (e < 0); the label c gives 
the transition. The inset shows a subset of the curves on a 
linear time axis; the increase of G(t) for short times cannot 
be resolved. 

Results of the (idealized) MCT equations for the po- 
tential part of the time-dependent shear modulus G(t) 
of hard spheres for various packing fractions <j) are shown 
in Fig. ^and calculated from fig : 

The normalized density fluctuation functions <&&(£) = 
{6g*(k,t)Sg(k))/ (\6g(k)\ 2 )/ enter, which are calculated 
self-consistently within MCT [25| . Packing fractions are 
measured in relative separations £ = (</> — 4> c )/4> c to 
the glass transition point, which for this model of hard 
spheres lies at 4> c — 0.516 ^5|. Note that this result de- 
pends on the static structure factor S(q) only, which is 
taken from Percus-Yevick theory, and that the experi- 
mentally determined value <f>'£ cpt ' = 0.58 lies somewhat 
higher [3, ■ For low packing fractions, or large negative 
separations, the modulus decays quickly on a time-scale 
set by the short-time diffusion of well separated particles. 
The strength of the modulus increases strongly at these 
low densities, and its behavior at short times presum- 
ably depends sensitively on the details of hydrodynamic 
and potential interactions; thus Fig. ^ is not continued 
to small times, where the employed model (taken from 
Ref. j24[) is too crude. For smaller (negative) separations 
from the critical density, e / 0, little changes in G(t) 



at short times, because the absolute change in density 
becomes small. Yet, at long times a process in G(t) be- 
comes progressively slower upon taking e to zero. It can 
be considered the MCT analog of the phenomenological 
Maxwell-process. MCT finds that it depends on the equi- 
librium structural correlations only, while HI and other 
short time effects only shift its overall time scale. Impor- 
tantly, this overall time scale applies to the slow process 
in coherent and incoherent density fluctuations as well as 
in the stress fluctuations |2f| . This holds even though e.g. 
HI are known to affect short time diffusion coefficients 
and high frequency viscosities differently. Upon crossing 
the glass transition bifurcation, the relaxation freezes out 
and the amplitude Goo of the Maxwell-process does not 
decay; the modulus for long times approaches the elastic 
constant of the glass G(t — > oo) — > Goo > 0. Entering 
deeper into the glassy phase the elastic constants increase 
quickly with packing fraction. 




FIG. 2: On the left side, increase of the Newtonian low 
frequency shear viscosity r/o relative to the high frequency 
one, 7700, for hard spheres as function of the packing frac- 
tion (f>; numbers (negative as cf> < <f> c ) denote the separation 
e = ((f) — <pc)/<f>c = ±W~^ /3 . Neglecting hydrodynamic 
interactions, ri-o — , , fcj E, T is obtained from the short time 
diffusion coefficient. The dashed line gives the asymptotic 
power law (4> c — (f)) 1 predicted by MCT with 7 = 2.46, while 
the relaxation times r, ; obtained from G(t^) = ksT / (2a) 3 
are multiplied by a factor 0.2 -7 and shown by a dot-dashed 
line. The inset presents the same data in a rectification plot, 
{Vo/v<x>)~ 1 ^~ 1 , where the power-law becomes a straight line. 
On the right side, the long time limits of the shear modulus, 
G(t — » 00) = Goo in units of ksT/a 3 , are shown versus <fi, 
with positive numbers n denoting the separation from 4> c . 

Calculating the viscosity according to Eq. @, i]o = 
Voo + Jo dt G(t) it is important to realize that the inte- 
gration is dominated by the slow process. The inset of 
Fig. ^ highlights this by showing G(t) on a linear axis 
of time. The viscosity increases as the area under G(t), 
which is almost exclusively set by the slowest process, 
which is the Maxwell one. The rapid decrease of G(t) to 
the plateau value Goo does not affect the viscosity as long 
as time-scale separation holds and the Maxwell process 
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is much slower. Figure |2] exhibits the obtained increase 
of the low frequency viscosity t]q for a quiescent colloidal 
dispersion of hard spheres. Also shown is the asymptotic 
power law predicted by MCT (with prefactor fitted), and 
the variation of the longest relaxation time t v of G(t) es- 
timated by G(t v ) — ksT / (2a) 3 . Asymptotically all three 
quantities exhibit the same dependence on </>, and diverge 
at the glass transition density <p c . The asymptotic power 
law variation holds for relative separations up to around 
-e «10 %. 
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FIG. 3: Storage part of the shear modulus G'(u>) correspond- 
ing to Fig. The inset shows storage and loss moduli (only 
for fluid states) for a number of densities. 
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FIG. 4: Loss part of the shear modulus G"(uj) correspond- 
ing to Fig. The inset shows the same data in a double 
logarithmic representation. 



easily. Storage G'(ui) and loss G"(uo) modulus are shown 
as functions of frequency in Figs. and 01 respectively. 
The slow Maxwell-process appears as a shoulder in G' 
which extends down to lower and lower frequencies when 
approaching glassy arrest, and reaches to zero frequency 
in the glass, G'(uj = 0) = Goo. The slow process shows 
up as a peak in G" which in parallel motion (see in- 
set of Fig. |3J) shifts to lower frequencies when e / 0. 
Including, in the spirit of the discussion in Sect. II, hy- 
drodynamic interactions into the calculation would af- 
fect the frequency dependent moduli at frequencies sat- 
isfying wthi ~ 1. For smaller frequencies, only a small 
correction would arise, G'^^uj <C 1/thi) ~ lutJoq and 
G' m (uj <C 1/thi) ~ (ujrjoo) 2 /G™, as long as time scale 
separation and consequently r/o 3> Voo holds. 

IV. SCHEMATIC MODELS FOR THE 
NON-LINEAR RHEOLOGY 

While within the MCT the linear viscoelastic regime 
has been worked out quantitatively, the recent general- 
ization of MCT to the non-linear regime 0, ^3] has not 
yet been fully solved. Only a number of universal pre- 
dictions have been obtained, which exist in any model 
that exhibits the generic bifurcation scenario from yield- 
ing solid to shear thinning fluid. In this section we verify 
that this scenario recovers the phenomenological Maxwell 
picture discussed in Sect. II. Hereto we use the most sim- 
ple schematic model that has also recently been success- 
fully employed to analyse flow curves from simulations 
0,|2i| and experiments [To|. 
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MCT thus provides a microscopic description of the 
linear viscoelasticity that follows the phenomenology first 
presented by Maxwell; of course in detail important and 
interesting differences exist, e.g. the exponential shape 
of the final decay postulated by Maxwell is replaced by 
a stretched-exponential decay in MCT. Tests of these 
specific predictions have greatly contributed to the veri- 
fication of MCT [3. 

While G(t) can be accessed by theory more directly, ex- 
perimentally its Fourier-transforms can be studied more 



FIG. 5: Normalized frozen-in glass parameters for a micro- 
scopic model, viz the elastic constant of a hard sphere glass 
Goo/G%>, and a schematic model, viz. the non-ergodicity pa- 
rameter f/fc of the Fi2-model, as functions of the separation 
parameter e; the normalization is the respective value at the 
transition. The upper axis gives e expressed in packing frac- 
tions for hard spheres. 

The central feature of the equations of motion in Refs. 
[Til fil| is that they contain the competition of two ef- 
fects. First, a non-linear memory effect increases with 
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increasing particle interactions ('collisions' or 'cage ef- 
fect') which leads to a non-ergodicity transition in the 
absence of shear, and second, memory effects vanish with 
time because of shear-induced decorrelation. Both these 
effects can be captured in the simpler 'schematic' models 
also. Importantly, the models can be made to obey sim- 
ilar stability equations as the microscopic approach, and 
thus the universal predictions hold. 

The well studied and comparatively simple schematic 
F^-model considers one normalized correlator $(i), 
which obeys a generalized relaxation equation: 

$(f) + T + jdt' m{t - t') = . (9) 

Without memory effects, m = 0, the correlator relaxes 
exponentially, <I>(i) — cxp —Tt, but with m/0, retarda- 
tion effects set in after a short-time variation (still given 
by the initial decay rate T, viz. <£>(£ — s- 0) = 1 — Tt + . . .). 
The correlator $(£) may be thought to correspond to 
the normalized shear modulus of the generalized Maxwell 
model of section II. A low order polynomial ansatz for m 
suffices to model the feedback mechanism of the cage- 
effect. We choose 

m(t) = i + JW 2 ( v±m + V2 * 2{t) } • (10) 

Without shear, this model has been studied extensively 
[THll2^ |. Increasing particle caging is modeled by increas- 
ing coupling parameters v%, vi > 0, and the only effect 
of shearing is to cause a time dependent decay of the 
friction kernel m. The system loses memory because of 
shearing. The role of the transport coefficient (viscos- 
ity) rj is played by the average relaxation time obtained 
from integrating the correlator, and this also is taken to 
determine the stress: 

/>oo 

a = 7 ? ; = 7 (r) = 7 / dt $(f) . (11) 
Jo 

At high shear rates, the memory function is strongly su- 
pressed, so that <J> returns to a single exponential, and 
the high shear visosity of the model follows as rj^ = 1/T. 

For the parameters of the model, yalues studied in the 
literature shall be taken [T3 . ITU [27j ■ The parameter T 
sets the time scale and determines the short time dy- 
namics. The bare Peclet number for the model thus is 
given by Peo = 7/r. Following earlier studies, the two 
interaction parameters are chosen as v-i = = 2 and 
vi = vf+s/l-y/v^-l), where tjJ = v§(y/4/v%-l) 0.828. 
Thus, a dependence on temperature or packing fraction 
enters the model only via s. A glass transition singularity 
lies at e = 0, where the long time limit <E>(i — > 00) = / 
jumps discontinuously from zero for s < to a finite 
value / > f c = 1 - 1/^/uf for e > and 7 = 0. The 
parameter / plays the role of the clastic constant G^ in 
this model. 

Comparing the amplitude / of the yielding process in 
the ivo ''-model with the elastic constants G^ calculated 
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FIG. 6: Correlators of the schematic -model as function 
of rescaled time tT. The curve marked with a relaxation time 
t taken at $(t) = 0.1 corresponds to a fluid state without 
shear, e = — ei = — 10~ 3,791 and 7 = 0. For increasing shear 
rates, log 10 (7/r) = -8, -6, -4, -2, 0, the correlators decay more 
quickly. The curve marked by the long time plateau value / 
corresponds to a glass state without shear, e = £2 = 16ei and 
7 = 0. For nonvanishing shear rates, log 10 (7/r) = -10, -8, 
-6, -4, -2, 0, the correlators decay more quickly to zero. The 
scaling function <j) + , whose integral gives the yield stress, is 
shown by a dashed curve shifted onto the si and 7/r = 10 -10 
curve. 

The inset shows flow curves of the model, viz. dimensionless 
stress a versus 7/r, for two states that would be fluid respec- 
tively glassy without shear; e = ±£2 and e = ±£3 = ±16£2- 
Straight lines with slope unity indicate the variation following 
the low shear 770 and high shear viscosity rjoa . 



for hard spheres in Sect. II. B, one notices a quantitative 
difference; see Fig. |5] While a change of e by 10 % causes 
an increase in by a factor 8, in the schematic model / 
changes only by 15 %. Fitting experimental flow curves of 
particles with strong repulsive interactions 01 thus the 
present F^— model can be expected to underestimate the 
increase of the stress values in the glassy region. 

The presence of a glassy arrested structure is equiv- 
alent to a frozen in part in the correlator or memory 
function; thus without shear $(t — > 00) = / > and 
m(t — ► 00) = g > hold for e > 0. With shear a non- 
decaying part in m(t) is impossible, as m(<7 3> 1) < 
(vi + V2) / ('jt) 2 ; as a consequence, also $(t) always de- 
cays to zero. Memory is cut off at long times, and Eq. 
(|10|l gives the most simple ansatz recovering this effect 
of shear advection in the microscopic equations, and the 
obviously required symmetry in 7. Figure shows for 
fluid and glassy, at 7 = 0, solutions of the F{ 2 -model 
the speed up of the relaxation caused by increasing shear 
rates. Integrating over the correlators, as given in Eq. 
pip, leads to the viscosity which consequently exhibits 
shear thinning; representative flow curves are shown in 
the inset of Fig. HJ 

It is interesting to analyse the origin of a yield stress 
more closely in the schematic model. This can be done 
most easily at the glass transition point, where the cor- 
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relator arrests at the (critical) value of plateau at long 
times, $(t — > oo ) = f c at £ = and 7 = 0; this quan- 
tity may be interpreted as the elastic constant GJ^of the 
amorphous solid at melting. Applying shear, and as- 
suming that the unique steady state has been reached 
after waiting sufficiently, the schematic model predicts 
a finite (critical) yield stress cr+ = dx&f (x) , that 
follows from the final decay of the (elastic) constant of 
the glass. The existence of a yield stress thus is seen to 
arise from the final relaxational process in the transient, 
which is only caused by shearing and whose relaxation 
time is given by 1/|7|. The phenomenology of the gener- 
alized Maxwell model of Sect. II. A is thus recovered in 
the schematic model, and can be recognized in Fig. El 



The part frozen in without shear of magitude /, decays 
during a process driven by and characterized by the same 
time scale as the external shear rate. 
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